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Abstract 

We introduce a class of distributions originating from an exponential family and 
having a property related to the strict stability property. Specifically, for two densities 
pg and p linked by the relation 

pg{x) = e'-c{e)p{x), 

we assume that their characteristic functions Jq and / satisfy 

fe{t) = r'^^\P{0)t) V t e R, 6* e [a,6] , with a and h s.t. a < < & . 

A characteristic function representation for this family is obtained and its properties 
are investigated. The proposed class relates to stable distributions and includes Inverse 
Gaussian distribution and Levy distribution as special cases. 

Due to its origin, the proposed distribution has a sufficient statistic. Besides, it 
combines stability property at lower scales with an exponential decay of the distri- 
bution's tail and has an additional flexibility due to the convenient parametrization. 
Apart from the basic model, certain generalizations are considered, including the one 
related to geometric stable distributions. 



Key Words: Natural exponential families, stability-under-addition, characteristic func- 
tions. 

1 Introduction 

Random variables with the property of stability-under-addition (usually called just stabil- 
ity) are of particular importance in applied probability theory and statistics. The signifi- 
cance of stable laws and related distributions is due to their major role in the Central limit 
problem and their link with applied stochastic models used in e.g. physics and financial 
mathematics. 
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In practical applications, stable laws appearing naturally as limit distributions for sums 
of random variables (the essence of the Central limit problem) are often used to describe 
the increments of stochastic processes. Stable distributions may exhibit a reliable fit of 
empirical data whose nature is related to random summation. However, certain peculiarities 
of the data (e.g., lighter tails than according to stable models) motivate the search for 
alternatives to stable laws. Among different ways to adjust the tail behavior of the stable 
distribution is the exponential (tilting), the smoothening of the tail of stable density is 
widely applied, see e.g. |14j . 

In the present work, we consider exponential families with a property that can be 
considered as a generalization of the usual strict stability-under-addition property (that is, 
with a real- valued function in place of the natural number of summands). The result is, as 
appears, related to exponentially tilted stable distributions. 

As the class we consider originates from natural exponential families and related to 
stable distributions, it is not surprising that our class appears to be a subclass of Tweedie 
distributions, that, for some combination of parameters, can also be characterized as ex- 
ponentially tilted stable distributions. 

Therefore, the result of our findings is a convenient representation of a certain subclass 
of Tweedie distributions - a characterization in terms of extended stability property. Like 
it is for the Tweedie distributions, the important advantage of the proposed class is the 
existence of a sufficient statistics which is due to its belonging to natural exponential 
families, - and the combination of the generalized stability with exponential tails. 



2 Definition of the class, derivation and properties of its chf 

Suppose X is a real random variable with cdf P and pdf p. A class of distributions associated 
with P is the natural exponential family with the pdf 

pe{x) = e'^c{e)p{x), (1) 

with the real parameter 9 G [a, 6], where a and b are such that G [o, With fe{t) 
denoting the characteristic function (chf) of pe{x) and f{t) the chf of p{x), the following 
equality is valid according to ([1]): 

/^(*) = 7P^' (2) 

The problem we investigate in relation to the above chf is the following: Do real valued 
functions a{9) and /3{6) exist, such that 

f9{t) = r^'\pm VtGM, ee[a,b]. (3) 

It is known (see e.g. |11| ) that if a chf is analytic in the neighborhood of the origin, 
then it is also analytic in a horizontal strip (either this strip is the whole plane or it has 
one or two horizontal boundary lines). It follows from the above that V G [O;^]) chf fg{t) 
is analytic {t G C) in the horizontal strip | Im (t) | < p{9) for some p{9) > 0. Then with 
the notation g{t) := ln(/(t)), one gets from ([2]) and Q the following balance equation 

git-i9) = g{-i9) + a{9)g{mt)- (4) 
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The above functional equation can be solved using an argument based on sequential dif- 
ferentiating w.r.t. to 6 and to t, which we postpone to Appendix. The solution w.r.t. the 
function g{t) has the form 

g{t)=A[l-{l-itcy], (5) 

where 7 < 2, A, c G M and the sign of A depends on the value of 7 (as we shall see, A < 
when 7 < and when 7 = 2; j4 > when 7 G (0, 1) and when 7 G (1, 2)). The expression 
for the chf f{t) immediately follows 

f{t) = exp{A[l-{l-itcr]}. 

Recalling ([2]), we can now write down the chf fg 



feit) 
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exp 



exp 



^•(1- 
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itc 
l-cO 



where B := 1 — c6. 

Note that the functions a{0) and P{6) are then explicitly expressed 

a{e) = (1 - c9y , (3(6) = 1/(1 -ce). 



(6) 



(7) 



2.1 Ranges of the shape exponent 7; role of parameters A and c 

The properties of introduced distributions differ with respect to various combinations of 
their parameters 7, A, c. 

As the role of the shape exponent 7 is crucial to relations between other parameters, in 
the present paragraph we analyze different cases that arise with respect to certain ranges 
of 7. That classification w.r.t. 7 also makes it easier to clarify the relation of the introduced 
class to stable laws and other important distributions. 

Returning to the representation for the chf / 

f^{t)=exp{A{[l-{l-itcy]} , (8) 

note that the properties of the corresponding distribution differ depending on the range 
of 7 (to mark the importance of this dependency, we shall use the notation instead of 
just /, and will point out the range of 7, such as f-y<:Q or /^g(;^ 2))- Listing different cases 
w.r.t. 7 below, we shall see that for each of the ranges of 7, is a proper chf only for a 
particular sign of A. In all cases, there exists a corresponding analytic function /^(t), t £ C, 
that agrees with the characteristic function /^(t), t G M. The domain of regularity for 
corresponding analytic function varies depending on the ranges of 7 considered below. 



Case a : 7 < Denoting 7 = — 7, for this case we have 

^<o(t) = exp {A [1 - (1 - iter'] } . (9) 

As 7(t) := (1 — itc)~'^ is the chf of a gamma distribution, we see that /^<o corresponds to 
the chf of the compound Poisson r.v. with gamma-distributed summands, provided that 
A:=-A>0: _ 

^<o(t)=exp{^[7(t)-l]}. (10) 
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Nit) 

In other words, /'y<o(i) is a chf of the random sum 5* = ^ Xi, where ~ Poisson and 

i=l 

Xi ~ gamma{'J, c) , i.e. 7 and c are the shape and the scale parameter correspondingly. 

Clearly, /'y<o(i)j t G C, is analytic when Im(t) G ( — 1/c, 00), with the entire line t = 
—i/c being the region of singularity for f{t). 

Case b : 7 £ (0, 1) This case, after a complex shift of the variable, corresponds to a 
certain subclass of stable distributions. Specifically, recall that if X = X + xq is a shifted 
one-sided Levy-stable variable then its chf can be represented as 

(j){z) = exp{izxQ — {—iz)°'a} , a > 0. (11) 

Clearly, (fTTj) with xq = brings us close to our case ([8]), differing by just a proper change 
of the variable. Specifically, with cz := ct + i we have 

/76(0,i)(i) = e^exp{-(l-itc)TA} ^ /^^(^^^^(^ _ i/c) = exp {-(-zzc)^A} , (12) 

which is close to (jlip up to the normalization factor (provided that A> Q). 

Furthermore, note that the re-parametrization cz = ct + i relates to the transformation 
of the density, similar to the one used in the introduction of /g, according to ([T|) and ([2]). 
When t = z — i/c then in order for /'y6(o,i)(^) to correspond to a proper chf we need 

a normalization factor /'ye(o,i)(~Vc)- From (|12p . we see that '/'^^"^^^^^(It/c')^ ~- f^^^'^i^) 
corresponds to one-sided stable chf and that /'yg(o,i)(~Vc) = e^. 

Hence, one-sided stable density corresponding to chf /^*"^ is the result of an exponential 
transformation (take the exponent parameter 1/c in ([1]) instead of 0) of the density corre- 
sponding to chf /-),g(o,i)- That is, the r.v. with chf /^g(o,i) i^^ turn relates to the exponential 
transformation of the one-sided stable density with the exponent parameter being —1/c. 
Note that both f^^""^ and /^g (0,1) are proper characteristic functions only when ^ > 0. 

Like in Case a, the function /'),e(o,i)(i) , t G C , is analytic in Im(i) G (—1/c, 00), but 
contrary to the previous case, there exists a limit for /-yg(o,i)(^) w.r.t t approaching the line 
t = —i/c (with the only point of non-regularity being (0, —i/c)). This limit corresponds to 
jstab introduced above which is the chf of one-sided stable law. 

Note that the above is in correspondence with one of the classical results on char- 
acteristic functions that says that a necessary condition for a function analytic in some 
neighborhood of the origin to be a characteristic function is that in either half-plane the 
singularity nearest to the real axis is located on the imaginary axis (see e.g. |llj for details). 

Remark 1 Check that the well-known Inverse Gaussian distribution as a special case with 
7 = 1/2. 

Case 0:7 = 1 This degenerate case corresponds to the constant r.v. whose cdf is the 
Heaviside step function, since the distribution corresponding to the chf e*"^'^* has a single 
unit jump at point Ac. 

Case d : 7 G (1, 2) Consider a re-parametrization 7 = 7/2 allowing to write down the 
chf f{t) as 

776(1,2) (0 = exp [a[1 - ((1 - itc)''f']] = exp [a[1 - (1 - cF)^]} = /~g(o,i)(-it), (13) 
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where t := 2it + ct^ . 

While /7g(o,i)(i) is a chf of the type considered above (Case b), /'yg(i,2)(^) is also a 
chf. Indeed, while </'(^,(j)(i) = exp — is the chf of the Gaussian r.v. (with 

expectation /x and variance cr^), the function 

/76(o,i)(-^*) = j e^''p^{x)dx = j [cl){2,2c){t)Y P^i.x)dx (14) 

which is a chf of some probabilistic distribution, as a continuous mixture of the Gaussian 
chf 4>{2,2c){t) = exp|2it — ct^} (provided that c := — c > 0) with the pdf p-^ whose chf 
/leCo.i) is the one from the above Case b. 

Since we already know from Case b that it should be that ^ > for /^g(o,i) to be a 
chf, in the present case we also have ^ > as a necessary condition for /^g(i^2) given by 
(|14p to be a chf. Like in Case b, /(t), t G C, is analytic only in the strip restricted by 
the horizontal line t = —i/c, i.e. when Im(t) G (— 1/c, oo), but additionally we have the 
condition that c < 0. 

Case e : 7 = 2 Check that /^=2 is the chf of the normal r.v. 

7^=2 (t) = exp {A {-2cit - ch^) } = 0(„2Ac,2Ac2) it) , (15) 
where A = —A. As the variance is cr^ = 2Ac^, it should be that A < 0. 



The sign of c Discussing the cases with respect to the parameter 7, we noted that the 
sign of the parameter A is crucial in each of the cases. Commenting on the sign of the other 
parameter, c, we noted that its sign is only crucial in Case d. In other cases, c can take any 
sign when is a proper chf due to the property of characteristic functions providing that 
/(— t) is a chf as soon as f(t) is a chf. However, we should keep in mind that the function 
f'y{t) (t G C), in all relevant cases, is analytic only in a strip depending on c. 

Limit behavior of / regarding mutual limit w.r.t. c and A Note that the one- 
sided stable distributions appearing in Case b and obtained via the change of the complex 
variable cz := ct + i could also be viewed as a limiting case w.r.t c — t- 00 and A ^ 0. 
Additionally, c and A should be linked so that A ~ c/c' (c G M), hence the limit w.r.t 
^ — > and c — )• 00 (keeping c constant) leads to 

Ae(o,i)(i) = e^exp |- Q - ztj (c^A)^ -— > exp {-(-it)^^} =: 4'^7oli)(*)' 

which corresponds to one-sided stable distributions appearing in Case b above. Recalling 
the notation of Case b, we see that /('*™)(t) = /^*"''(i) = exp{-{-itc)'^A}. 

2.2 The role of the natural exponent 9 

Let us turn again to the representation for fg. Rewrite ^ noting that an additional re- 
parametrization with A := A ■ and cg = leads to the form ([8]) 

f^,e{t) = exp {I - [1 - (1 - itcgy]'} , (16) 
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where we used the notation f^^g instead of just fo, indicating the significance of both of 
the parameters. 

Clearly, Cases considered above in relation to chf can be revisited for fg, but apart 
from the ranges of the shape exponent 7, the role the natural exponent 6 is important. 

Obviously, the borders of the analyticity strips depending on c, as well as the conditions 
on A, will change in correspondence with 0, i.e. one can view introducing f^^g in place of 
as passing from 9 = to non-zero 9, so that the "new" c and A will depend on 9. 

In particular, replacements through Case a - Case e are the following: 

• the strips of analyticity: (— i/c, 00) =^ (— i(l/c — 9), 00) in Cases a, b and d; 

• the signs of coefficients: recall that Case a and Case e imply A < 0, while in both Case 
b and Case d we have ^ > 0; for /^g the same holds with A =^ A-B'^ = A-{l — c9)'^. 

Revisit, for instance. Case d where we have a representation analogous to the one 
expressed by (fT3]) and (fn|) : 

I-fe{i,2),e{^) = fie{OA),e (2* " ^^i^) = J [^i2,2ce)i*)VPiA^)dx, (17) 

where pdf p^ g corresponds to the chf /i^eCO,!),^ = ^xp ^A [l — (1 — itcg)^'^ | and (p{2,2ce) 
is a chf of the normal r.v., i.e. 4>{2,2cg){t) = exp |2it — cgt^} (with cg := —cg). So that in 
order for /^g(i,2),6' to be a chf, it should be that cg = c/(l — c9) < and A = A ■ B'^ > 0. 

All other cases w.r.t 7 can be revisited for fry^g analogously. Additionally, certain values 
and ranges of 9 correspond to particular special cases which we consider below. 

Range of 9, singularities and corresponding limits. While c can be positive or 
negative, the point when ^ = 1/c is the case of singularity and/or the limiting case. 

Specifically, Case a with 9 = 1/c is the case of singularity of the chf f^^g. In Case b, 

the limit 1/c ^ 9 exists with limiting chf's f^^'^\t) = exp { — (— izc)'''^}. 

In a way, approaching the point 1 /c with respect to 9 means approaching correspond- 
ing stable (or mixed normal-stable) distributions, while increasing the exponent parameter 
9 corresponds to lightening the tails of distributions Pg to exponential ones. 

Limit case 9 — > 00 . The question in interest is: what distributions appear in the 
limit case 9 — > 00? In |T|, limit laws for the whole class of natural exponential families 
w.r.t. the growing exponent parameter were investigated. It was shown that for the r.v. 
Xg with density of the form ([T]), constants ag > and hg can chosen so that a limit 
{Xg — bg)/ag — > Y, 9 ^ oo exists and that only possible limit distributions Y are the 
normal and gamma distributions. Specifically, the following results were proved. 

• If Pg, 9 £ Q, is the cdf of the exponential family Xg whose pdf is given by ([ip and if 

there exist constants ag > and fea € R such that Pg ( 2—^ ) 5. G weakly to some 

\ "■<> J g^^ 

non- degenerate cdf G then G belongs to the so called extended gamma family, which 
includes normal and gamma densities. 

• In the case of convergence to a non- degenerate limit G, ag and bg should be taken 
such that Pg is centered and scaled by expectation and standard deviation. 
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Criterion of the convergence of an exponential family ([T]) to gamma/normal distribution 
can be expressed in terms of the moment generating function (mgf) associated with p{x), 
i.e. M(A) = / e^^p{x)dx. Specifically, as proved in |Tj, the following holds. 

If the function s(A) = 1/V m"{X), where m(A) := InM(A), is self-neglecting, then the 
exponential family Xx, X £ Q, is asymptotically normal. Otherwise, the only possible limit 
distribution is gamma distribution. 

It is easy to check when the normal distribution appears as a limit in our case. Recalling 
([5]), for the log-mgf (or cumulant generating function) of the introduced family we have 

m(A) = A[l-(1-Ac)^], 

so that 

m"(A) = ^0^7(7 - 1) [(1 - Ac)^^2j ^ 

Clearly, the above condition on s = 1/V m" is only satisfied when 7 > 2. That means that 
in all the cases w.r.t. to 7 that are relevant for our study, the only possible limit distribution 
G is gamma distribution, i.e. according to the usual terminology, our distribution lies in 
the domain of attraction (DA) of gamma distribution. 



3 Relation to the Tweedie class and special cases 

Like briefly noted above, the proposed distributional family includes some well-known mod- 
els, which can be considered as special cases of exponential family with stability property. 
At the same time, it itself is subclass of another class ~ the Tweedie distributions [?], which, 
in its turn, is a part of even a wider class of exponential dispersion models. 



3.1 Relation to the Tweedie class 

For general definitions and characterizations of exponential dispersion models we refer to 
[?]. The pdf of an exponential dispersion model (EDM) is represented in a way similar to 
that of our class 

fEoiy) = Pi3{y) exp i^^{9y - K{e)) 

where /3 > 0. A convenient representation of the Tweedie class is through the moment 
generating function M{t) = f exp{ty)f{y)dy. The cumulant generating function of the 
EDM is 

TTlEDit) = log MEoit) = [K{e + tc/)) - k(0)]//3 . 

As the derivatives of the function k(-) w.r.t 9 provide the values of the cumulants of the 
distribution, k(-) is called the cumulant function. So that we have fj, = k'{9) for the mean 
and (3k," (9) for the variance. Since the mapping from 9 to fi is invertible, there exists a 
function V{-), called the variance function, such that k"{9) = V{fi). 

The Tweedie class is then defined as EDM with the variance function V{-) of a specific 
form, namely V{ii) = ijP for p G (— oo,0] U [l,oo) and V{fi) = exp(t;//) for p = 00 (here, 
V ^ 0). The correspondent cumulant function k for the Tweedie class can be represented 
as 



when p ^ 2 



2-P ' ^"""'^ y ^ ^ (18) 
log , when p = 2 . 
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Defining a new parameter a via 



{p-m-a) = 1, 

a representation for tlie moment generating function of tlie Tweedie class is obtained: 



' exp{i^=gf [(l + s/^r-l]} , when p / 1,2, 

(1 + sBy^lf^ , wlien p = 2, (19) 
exp |-|e^(e* — 1)| , when p = 1. 



From above, comparing with the chf ([5]) of our class, it is clear that our class is actually 
a subclass of Tweedie distributions with some parameter p ^ 1,2 and a certain combination 
of other parameters. 

Correspondingly, the basic properties that we discussed above do hold for the Tweedie 
class as well. For instance, it is known that the Tweedie distributions can be obtained 
through exponential damping of stable distributions. 

The distributions of the Tweedie class are known to have a scale invariance property, 
that can be expressed through the pdf fj-w as 

fTw{x; /i, (p) = cfTw{cx; cfi, c^'^fS) , for c> , 

which also resembles stability property (in terms of scale invariance), however is not directly 
related to the generalized stability property introduced by considering the relation ([3]). 

Note also that our subclass does not explicitly contain the parameter /3 in the above 
ED notation (parameter related to the variance of the distribution). 

Besides, note that gamma distribution belonging to the Tweedie class, in our case 
appears as a limit case — t- oo, like mentioned in the previous section. 

3.2 Inverse Gaussian and Levy distributions 

Two important special cases are visible immediately from the form of the characteristic 
functions and fry g. 

Recall the characteristic function representation of the well-known inverse Gaussian 
distribution 



9lG{t) = e" 



1/2 



(20) 



and note (as already briefly remarked in Case b of Paragraph 12. ip that it can be viewed as a 
special case of the introduced exponential family with stability — just plug the parameters' 
values {A = = ^a~'^ ~ iii^° general representation ([S]). 
Recall an important property of the inverse Gaussian distribution: 

• If Xi has an IG{fiWi, Xwf) distribution for i = l,2,...,n and X/ are indepen- 

n 

dent, then S = Yl'i=i-^i ~ IG{fj/w, Xw'^), where w = '^i- other words, X ~ 

i=l 

IG{fi,X) ^ tX ^ IG{tfi,tX) 
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Note that this property is related to the scale invariance property and recall that stable 
distributions are also included in certain subclasses of introduced exponential family (see 
Case b and Case d), so that the scale invariance appears to be naturally embedded in our 
model. 

Another well known distribution appears as a limiting case of the Inverse Gaussian 
distribution: the chf fievy of the Levy distribution is clearly a limit of gjc w.r.t. fi — > oo 
(and respectively A — > 0) which leads from (j20p to 

fLevy{t)=e-^^'. (21) 

Inverse Gaussian and Levy distributions are the important cases as they have simple ex- 
plicit forms for their pdf's. 



3.3 Note on Tempered Levy processes 

The combination of the properties of our subclass relates in its purpose to other modified 
models such as the so-called Tempered Levy process. 

The idea behind Tempered Levy processes (also known as Truncated Levy flights |14| ) 
is to use a distribution which coincides with the stable one around zero (i.e. for small 
fluctuations) and has heavy tails yet decreasing to zero fast enough to assure a finite 
variance. 

The latter is achieved by truncating stable distributions, and one of the ways the trun- 
cation can be done is the so called smooth exponential truncation. The result corresponds 
to an infinitely divisible distribution enabling an analytical chf representation: 



(t) = exp 



cos(7ra/2) 



(fi^^ + t^)°/^ cos ( a arctan ( — — /u" 



where is the cut-off parameter. The asymptotic behavior of (p is given by 

(pa; flit) ^ ~_ (j)a(t) 



where (pa{t) is the chf of the (non-truncated) stable distribution, i.e. for small values of x, 
the pdf pa-fj_{x) corresponding to the chf 4'a;fi behaves like a stable law of index a. 

This concept is then interpreted in terms of Levy processes and associated Levy mea- 
sures. Specifically, if i^(x) is a measure of a Levy process (generally, defined as the expected 
number, per unit time, of jumps whose size is less than x), then the measure 

u{dx) = e~^''u{dx), 

corresponds to another Levy process, whose large sizes are "tempered" with exponential 
damping. 

The ideas are therefore similar to the one that lead to the introduction of the Tweedie 
class, as well as to the characterization of our subclass, yet with certain methodological 
difference: The idea of exponential tempering utilizes stable distribution as the initial dis- 
tribution whose tail is then exponentially smoothened, while in our approach the stability 
property appears naturally from the properties of the introduced class. That allows to 
view the exponential family with stability property as a natural extension of both classes: 
stable distributions and exponential families which is promising regarding their practi- 
cal applications, including the use of sufficient statistic naturally inherent in exponential 
families. 
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4 Modifications related to geometric stable laws 



Preliminaries and definitions First, recall the origin of geometric infinitely divisible 
(GID) distributions. 

The chf of an infinitely divisible (ID) r.v. X implies the representation (f>{t) = ((/>„(t))"', 
where the chf (f>n{t) corresponds to some r.v. X^'^\ which in turn implies that the sum 
X[-^ + ... + xt^ of iid copies of has the same distribution as X. The so-called 

transfer theorems (see e.g. [3|) state that the random sums 

+ + (22) 

where (i/„,) is a sequence of integer- valued r.v.'s such that f„ — ^ oo (in probability) while 
i^n/n- (in distribution), converge (in distribution) to a r.v. Y whose chf has the form 

u}{t) = C{- In (pit)), (23) 

where C is the Laplace transform of the r.v. v. 

If Un = fp has a geometric distribution with mean 1/p (where p is close to zero, and 
i^p is independent of X'^s) which converges in distribution to the standard exponential 
distribution with Laplace transform C{u) = (1 + u)~^ then (|23p turns into 

io{t) = {l- In cl){t))~\ (24) 

and random variables whose chf has the form were introduced in [6] as geometric 
infinitely divisible (GID) random variables. Since sums such as (|22p frequently appear in 
many applied problems, GID have a variety of applications. The geometric stable distribu- 
tions (GS), also originated from [6] and developed in later works e.g. [3] and [5], appear as 
a natural subclass of GID and are widely used for modeling of stochastic processes. The 
chf of GS laws has the form (|24|) with (j){t) being the chf of a stable r.v. 

It is natural to extend the exponential family with stability property introduced in 
previous sections w.r.t geometric stability. Specifically, for the chf g of the form (fT6|) . 
the corresponding geometric extension will have the chf 

'^lAt) = 1 1 \ V t G K, ee[a,b]. (25) 

Since /■y^g{t) = exp |y4 • [1 — (1 — itc5i)'^]|, note that when A = A • B"' < (which corre- 
sponds to certain ranges of 7) then (|25|) can viewed as a sum of geometric progression with 
the initial value (1 — A)^^ and common ratio (1 — itcg)^ . 

Revisiting the cases considered in Paragraph 12.11 again, this time in relation to uj-yfi, 
the argument above leads to a convenient interpretation of Case a, which we sketch below, 
along with other cases w.r.t the ranges of 7. 

• Case a 7 < 0. In this case, /-y<o,6» is a chf of a probability distribution only when 
A <Q, which coincides with the condition under which the representation (|25p can be 
interpreted as a geometric progression. Thus the cdf VL^^q q corresponding to the chf 
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<^7<0,6 the representation as a series of convolutions of the gamma distribution 
with itself: 

n,^,^,{x) = -—jY. (26) 



n=0 



where F^fi{x) = 1 — F^fi{x)^ and Fiyfi is the cdf of gamma distribution with shape 
parameter 7 = —7 and scale parameter cq. 

Case b 7 G (0, 1). In this case /'ye(o,i),6' is a chf only when ^ > 0, so the representation 
of the form (|26p is not valid. However, as the representation (j25p implies that any 
oj^fi is a chf as soon as is a chf, it follows that ti;^g(o,i),0 is the chf (provided 
that ^ > 0), and the corresponding r.v. is a geometric analogue of exponentially 
transformed one-sided stable r.v., corresponding to the chf /^g(o,i),e- 

Case d 7 G (1,2). Recall the representation (p7|) meaning that /'yg(i,2),0 is a chf 
linked via the transformation j \_4'{2,2ce){^)Y^ Pifii^)'^^ with the pdf p^y^g of Case b 
(7 G (0, 1)) and with the chf 4>(^2,2cg) of normal distribution (with mean and 
variance (2,2cg)). Then a;^g(i,2),6» given by (|25p is also a chf, the geometric analogue 
of /'ye(i,2),0i provided that A > and additionally cg < 0. 

Case e 7 = 2. As chf /■y=2,e corresponds to the Normal r.v., uj^=2,d is the chf of the 
geometric analogue of the Normal r.v. 

Check that this case can also be interpreted in the way similar to the Case b, so that 
10^=2^ can be viewed as an exponentially transformed Laplace r.v. 

Case f 7 > 2. This is a new case in our framework, since /^>2,6» (t) does not corre- 
spond to a proper probability distribution (in the representation J [e^*-*^] pg{x)dx 
that we used in Case b and Case d, the polynomial P{t) can not have a degree 
greater than 2, due to Marcinkiewicz's theorem |12j). 

However, a;^>2,e(i) = (1 — In /7>2,e(i))~^ ^-^ ^ proper chf. That could be seen if we 
assume 7 G (2,4), denote 7 := 7/2 and a := 1/(1 — A) and consider 
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^76(2,4), 6» (t) 



l-A{l-itce)^ l + A[{l-itcgp]^ 1 - {I - a) [i-izcep]'^ 

with z = t+i/c. Clearly, the r.h.s. of the above can be represented as Cl{— log f~^°'^{z)) 
where f~'°'^{z) = exp{y4[— (— izc^)'^]} is the chf of one-sided stable random vari- 
able, as Cl{u) = 1/(1 — b'^u'^) is the Laplace transform of the Laplace distribu- 
tion. Clearly ujg{z) := 1/ ^1 — (1 — a) [(— izcg)^]^^ is a chf of geometric analogue of 
one-sided stable r.v. (provided that (a — 1) > 0), and since ujg{—i/c) = 1/a, then 
w^g(2,4),e(0 = ^g{z — i/c)/ujg{—i/c). Then a;^g(2,4),e is itself a chf. Note that the 
condition (1 — a) > means that it should be that A> 1. 

When 7 > 4, we can argue analogously. Note that the r.v. corresponding to the chf 
^■y>2,e is nevertheless not geometric stable, nor it is geometrically infinitely divisible. 
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5 Summary, interpretation and further extensions of the model 



5.1 Summary of properties 

As the classification of different cases witfi respect to tlie ranges of 7 spreads out to f^^g 
and uj.y^g, it might be useful to summarize the relevant information in a summary table. 





Case a 
7 < 


Case b 
7G (0,1) 


Case d 
7G(1,2) 


Case f 
7 > 2 




exp{:4[7(t)-l]} 
• A<Q 


with z := t + i/c: 
exp { —[—izc]"'' A } 
• A > 


with z := t + i/ c: 
exp{-[(-izc)2]"yl} 
• A > 






g-I[7(t)-l] 

• I< 


with z := t + i/c: 
exp 1 —[—izcQ\^ a| 
• I> 


with z := t + i/cg: 
exp{-[(-izce)2pl} 
• I > 




Geometr. 

^1,9 


00 
n=0 

• A < 


(l-ln^e(o,i),8W)"' 
• ^ > 


(l-ln/^6(i,2),eW)~' 
• A > 


1-6 

l-b[{-izce)'iY 
• I> 1 



5.2 Interpretation in terms of a subordinated Levy process 

While in Section H] a new class was introduced whose chf (p is linked with the chf g 
through (f){t) = C{— In fy^g{t)), the chf f^^g itself, for 1 < 7 < 2, can be represented in the 
above form. Indeed, according to Case d of Paragraph 12.11 we have the following mixture 
representation for /^g(i^2) 

/7e(i,2)(i) = j [(Pi2,2c)it)]'" P^&io,i)ix)dx = J e"'°<^(2.2.)Wp7e(o,i)(2;)c?2; = £^(-ln0(2,2c)(t)), 

(27) 

and the same is valid for /-yg(i,2),e ^ well. 

Recall that a r.v. 5 with chf of the form C]sf{— In fx{t)) can be interpreted as a random 
sum (with random but finite number of summands), if the Laplace transform Cm corre- 
sponds to a discrete random variable (a realization of a counting process). One of the most 
celebrated examples is the Poisson process whose Laplace transform is exp |A (e~* — l)}, 

N{t) 

SO that the chf exp {A {4>{t) — 1)} corresponds to the compound Poisson r.v. S = ^ Xi. 

i=l 

According to (|13p . in our case the LT of the r.v. associated with the counting process is 
^n{u) = exp 1^4 [1 — (1 + cu)'^~\ } with < 7 < 1. Clearly, it corresponds to a continuous 
r.v.; moreover, as discussed in Case b, the chf /::yg(o,i) of this r.v. is obtained from the one- 
sided stable chf by a complex shift, s.t. the corresponding r.v. relates to the exponential 

transformation of one-sided stable density. Hence the chf (j27p cannot be interpreted as a 

N{t) 

chf of a compound r.v. in its usual form ^ Xi. 

1=1 

However, the corresponding r.v. can be interpreted in terms of continuous-time ana- 
logues of random sums — the increments of a subordinated Levy process. For the chf of any 
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Levy process Xt on IR we have (see e.g. [2]) 

where ljJx^ usually called the characteristic exponent of the Levy process, is defined via 
4>x{u) = e'^-^^") , with (px being the chf of the r.v. Xi . 

Furthermore, for a subordinated Levy process Yf = with subordinator S{t) (which 

can be viewed is a random compression of time) we have 

E [e"^'] = e*- isi^xiz)) ^ ^ g ][j ^ (28) 

where ls{z) defined via E |^e"'^^*^] = e*''^^"^ is called the Laplace exponent of S. For the chf 
of the increment Yi in unit time, ()28p turns into e^^^'^^^^^\ 

Recall that our chf /'yg(i^2)) according to (f27|) . can be represented as 

f^eii,2){z) = ln</.(2,2c)(^)) = exp [a [l - (1 - ln</.(2,2c)(^))^] } = e'^0"'^(2^^-)(^)) , 

where It^ is the Laplace exponent of the Levy process whose increments have the chf /::j;g(o,i)i 
i.e. f^(^[Q^i){iu) = C^{—u) = 6'^*-"^ Therefore /^g(i^2) can be interpreted as a chf of the r.v. 
Yi, i.e. of the increments of subordinated Levy process Yf = ^5{t) iu unit time t = 1. The 
increments of the underlying process Xt have normal distribution with parameters (2,2c), 
while the increments of the subordinator S{t) are distributed according to the exponentially 
transformed one-sided Levy distribution whose chf is given by /;^g(o,i)- 

A possible way to extend the model is to use any relevant chf as the component (j) in 
the representation (|27p . That would also lead to an interpretation of the corresponding 
distribution as a distribution of the increments of a subordinated Levy process Yt = 
with the same subordinator S{t) but with a different underlying process Xf. 



A Derivation of the chf in the expHcit form 

The function / to be found satisfies the balance equations 

feit) = ^TT^ ^^^^ fsit) = r^'HPm V t G E, ee[a,b], Oe[a,b]. 
In terms of the log-function g{t) = log/(t), it turns into 

g{t-ie) = g{-ie) + a{9)g{mt)- 

In order to prove that g can be represented in the form ([5]), let us differentiate g{t — iO) 
w.r.t. to t and to 9. First differentiate it w.r.t. to t 

g'{t-ie) = a{e)-m-9{mt)- (29) 

Considering a limit w.r.t — )• and assuming that a(0)/3(0) ^ 0, we get 

g\t) = Ag'ibt), (30) 

where b := /3(0) and A := a(0)/3(0). Two major cases are possible: b ^ 1 and 6 = 0. Let us 
start with the first one. 
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6^1 

Without loss of generality, assume that \b\ < 1 . Let us differentiate (j30p further k 

IjllXlGS 

5(fc+i)(t) = A.6V=+i)(6t), 

and consider it at t = 

Considering a limit w.r.t k — t- oo, we see that 1 < Ab^ — t- 0, which means that 
^('=+i)^Qj = when k > Kq (for some large enough Kq) . 

It follows from the latter that g{t) is some polynomial g{t) = P{t), so that the chf /(i) is 
the exponent of a polynomial f{t) = e^^^\ Due to Marcinkiewicz's theorem |12) . we can 
conclude that the degree of P should not be higher than 2, e.g. it is a quadratic polynomial. 

In other words, h ^ I 'm. (j30p corresponds to the case when f[t) is the chf of the Normal 
distribution considered in Section 2.1 as Case e. 

6 = 1 

Specifically, 6=1 means that not just /3(0) = 1 but also a(0) = 1, as follows from ()29p . 
Denoting A(0) := a{0)P{6) and keeping in mind that A(0) = 1, let us now differentiate 
(|29p w.r.t to 6 which gives 



(31) 



-ig"{t - iO) = ^{e)g'{mt) + mP'{0)g"{mt)t 
and considering a limit w.r.t 6 — t- 0, we get 

-ig'\t)=h'{Q)g'{t) + tp'{Q)g"{t), 

which leads to 

-{tfi\Q)+i)g"{t)=A'{Q)g'{t) 
Denoting u{t) := g'{t), we get an equation 

u'{t) _ -A'jO) 

u{t) ~ {i + t/3'(0)) 

The solution of this equation gives 

Inu(t) = -ci ln(i + t/3'(0)) + C2 u(t) = _^ -^^^ (where c := e' 

i.e. g'{t) = c(t/3'(0) + which means that 

»<"= (l-c°),3-(0) W°)' + i) 

Denote 1 — ci =17 and, in order to express the constant C2 explicitly, consider the above 
at t = 

c . 

-l^ = -C2, 



C2^ 



(1 - ci)/3'(0) 

so that we can denote Ai := and rearrange to get 

g{t)=A, [{p'{0)t + ir-i^]. 
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With two more notations A := —AiP and c := /3'(0), after rearrangement it finally turns 
into 



Throughout the derivation, we assumed that /3'(0) 7^ 0. Let us now consider the opposite 
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g{t) = A[l - {1 - icty] . 



case. 




282-289. 
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